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We present recent results on the Landau gauge gluon and ghost propagators in SU (3) pure gauge 
theory at Wilson j3 = 5.7 for lattice sizes up to 80'* corresponding to physical volumes up to 
(13.2 fm)'*. In particular, we focus on finite-volume and Gribov-copy effects. We employ a 
gauge-fixing method that combines a simulated annealing algorithm with finalizing overrelax- 
ation. We find the gluon propagator for the largest volumes to become flat at 0.01 GeV^. 
Although not excluded by our data, there is still no clear indication of a gluon propagator tending 
towards zero in the zero-momentum limit. New data for the ghost propagator are reported, too. 
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1. Introduction 

Presently, there is an intensive exchange of results and opinions between groups carrying out 
analytical and numerical studies of infrared QCD. This ongoing research focuses in particular on 
the infrared behavior of the Landau-gauge gluon and ghost propagators. The latter is intimately 
related to the confinement [[I], |3|]. What makes analytical predictions possible also in the non- 
perturbative sector of the theory is the possibility to write down a hierarchy of Dyson-Schwinger 
equations (DSE) connecting propagators and vertices. Under rather mild assumptions the hier- 
archy can be truncated. However, not all assumptions have been thoroughly checked. Note that 
recently the full system of Landau gauge DSE has been solved without any truncations within the 
asymptotic infrared region with a power ansatz for all Green functions involved Numerically, 
the propagators can be studied from first principles in terms of Monte Carlo (MC) simulations of 
lattice QCD. It is worth to compare the lattice results with the asymptotic power-like behavior, and 
with numerical DSE solutions found in finite volumes [||]. It is interesting to see whether there 
remain differences as the infinite-volume limit is approached. 

On the lattice we approximate the gluon propagator as the MC average 

D%{q) = m)A'A-4)) = 5"' (5,. - ^) ^ (1.1) 

\ H J H 

with the gluon field i4^-^^/2,/j = (1 /2iago)(^7v,^ — ^7rV)tiaceiess transformed into Fourier space. The 
lattice momenta =2 n k^/L^ with integer k^ £ (— +L^/2] are related to their physical 
values by = [2 /a) sin (tt k^/L^). 

The ghost propagator in momentum space at non-zero q^ is defined by double Fourier trans- 
formation 

&^\q) = L-Mx-y) [M-'fb\ = §-bZghip_ ^j_2) 

x.y ^ ' ' ^ 

Practically, this is done by an inversion of the Faddeev-Popov (F-P) matrix M^^ using a conjugate- 
gradient algorithm with plane waves as sources. The Faddeev-Popov operator in terms of the 
Landau gauge-fixed links is 



(1.3) 



with T" = X" /2 {X" are the Gell-Mann matrices). The functions Zgi{q^) and Zgi,{q^) are called 
dressing functions of the respective propagator. 

In Landau gauge the gluon and ghost dressing functions are predicted to follow the simple 
power laws [|6|] 

ZgH{q')oc{q^)-'^ and Zg,{q^) o. [qY'^ . (1.4) 

in the asymptotic regime q^ 0. Thereby, both exponents are related to some K which, under 
the assumption that the ghost-gluon vertex is infrared-regular, takes a value of about K = 0.596 
[^]. That is, the gluon propagator is predicted to decrease towards lower momenta and to vanish 
at q^ = 0. At which scale this asymptotic behavior sets in cannot be concluded from those studies, 
however. 
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A truncated system of DSE formulated on a 4D torus and numerically solved [|5|] predicts 
a specific finite-volume behavior which, at the first glance, looks quite similar to earlier lattice 
results obtained in particular by some of us Characteristic deviations for momenta ^ ~ 1/L of 
the gluon and the ghost propagator from the momentum dependence of the respective propagators 
at infinite volume should be expected. In order to check the DSE predictions we decided to evaluate 
the gluon and the ghost propagator for increasingly large symmetric lattices. We have measured 
the gluon and ghost propagators for configurations generated with the Wilson gauge action at fixed 
j3 = 5.7 on 56^, 64^, 72^ and 80"^ lattices. Note that the latter corresponds to a volume of about 
(13.2 fm)^. Comparing results from either analytic or numerical approaches for varying 4- volume 
will hopefully allow (i) to conclude for which momenta data on both propagators are reliable, and 
(ii) to estimate the order of magnitude of distortion of the momentum dependence due to finite-size 
and Gribov-copy effects. This paper presents first results of this study. 

2. Gauge fixing 

To fix the Landau gauge, we apply to all links a gauge transformation g £ G (G = SU (3)) 
mapping Ux^n ^Ux,^ = gxUx.^gl_^_^, with the aim to maximize a gauge functional 

^f^[^] = T^I^^Tr^[/,,^, (2.1) 

or, more exactly, to find the global maximum of Fjjig] [|l|]. In this work gx £ G is considered as a 
periodic field on all lattice sites. To find the global maximum in practice is a complicated problem 
which becomes exceedingly time consuming with increasing lattice volume. Starting from an initial 
random gauge transformation gx one generally arrives at one of many local maxima of Fu [g] . The 
corresponding gauge-fixed configurations are called Gribov copies. They all satisfy the differential 
gauge condition d^A^ = together with the additional necessary condition that the Faddeev-Popov 
operator has a positive spectrum (apart from its 8 trivial zero modes). With increasing volume, the 



copies become dense with respect to the value of the functional ( |2.1| ) and the spectral density of the 
F-P operator near zero grows [^. 

One way to suppress the effect of the Gribov ambiguity is to find A^copy local maxima of Fu [g] 
and to choose among them the "best" one ("be"), which possesses the largest value of Fu[g]. The 
underlying idea is that the maximal value of the local maxima approaches the global maximum of 
F[/[g], and the distortion of gauge-dependent observables, computed on such copies, vanishes in 
the limit A'^copy °°- Such studies normally use the overrelaxation (OR) technique to search for 
the maximum of Fij[g]. They have been carried out in |j8|, |l^, 11, 12 1 and have shown the Gribov- 



copy effect to become weaker with growing lattice extension L - in accordance with Zwanziger's 
conjecture [P. This suggests that it is tolerable to restrict gauge-fixing computations on large 
lattices (L > 48) to one gauge copy only^ Our simulations of the SU (3) propagators at L = 56 
were carried out using an OR algorithm with the overrelaxation parameter set to a = 1 .70. The 
number of gauge-fixing (GF) iterations did not exceed 10"* in most of the cases. However, we 
find a considerable slowing down of the OR GF process on a 64^ lattice. This was one of the 
reasons why we switched from using OR to a simulated annealing (SA) algorithm. SA, also known 



This is called "first copy" ("fc") in a multi-copy approach. 
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as a "stochastic optimization method" has been proven to be highly effective in coming close to 
the global maximum in various problems of different nature with multiple local maxima. It was 
proposed in [13, 14] and has found numerous applications in various fields of science. The idea 
of applying SA for gauge fixing has been first put forward and realized in the case of maximally 
Abelian gauge in [15]. The method is designed to keep the system long enough pending in a region 
of simultaneous attraction by many local maxima during a quasi-equilibrium process undergone by 
the "spin system" formed by gx interacting through the fixed {f/r,^} field with Fu[g] as energy. The 
temperature T of the spin system is decreased by small T-steps between updates in a range of T 
where the penetrability of functional barriers strongly changes. Theoretically, when infinitely-slow 
cooling down to T = 0, the SA algorithm finds the global maximum with 100% probability. For 
complicated systems with large numbers of degrees of freedom and of functional local extrema, 
e.g., for GF on large lattices, we have to restrict the number of A^iter T-steps of cooUng the system 
from Tmax to Tmin to, Say, 0(10^^). Within these limits we can still try to attain an as high value of 
the functional studied (in our case, as possible. Note that GF with SA requires a finalizing 

OR in order to satisfy transversality d^A^ = with a given high precision. 

For the SU (3) case we chose T^ax such that it leads to a sufficiently large mobility in the 
functional space. The final temperature Tmin was taken low enough that the subsequent OR was not 
slowed down while penetrating further functional barriers. This is witnessed by the check that the 
violation of the differential gauge condition, {d^A^Y, monotonously decreases until the machine 
precision is reached (stopping criterion) in almost all cases. In practice, for L = 56,64,72 and 80 
we restricted ourselves to one copy, and carried out from 5 x 10^ to 15 x 10^ heatbath (HB) sweeps 
of SA with 4 microcanonical sweeps after each HB one. We checked that the smaller T-steps 
are done in between, the higher the local maxima being reached. Finally, we note that a linear 
decrease in T seems not to be the optimal choice. T-schedules with smaller T-steps close to T^ax 
and larger T-steps at the end (with A^iter fixed) lead to higher 7^[/[g] -values (after completing the full 
SA procedure). 

3. Ghost propagator results 

The SU{3) ghost propagator at j8 = 5.7 is shown as a function of in Fig. [l]; on the left 
hand side for a single 56^^ configuration, simply comparing results after either OR or SA gauge 
fixing, and on the right hand side as an average over 14 configurations in the case of OR, and over 
7 configurations using SA gauge fixing. The influence of the gauge-fixing method, here through 
the emerging copy, can be seen only for the three lowest momenta at this lattice size. In general, 
one notices that the higher the gauge functional, the lower the estimates of the ghost propagator at 
the smallest momenta. This comparison (Fig. [l]) demonstrates that the problem of Gribov copies 
does still exist for L = 56, resulting in maximally 10% difference of ghost propagators at lowest 



momenta. Note that this result cannot be directly compared to our previous studies [ ]10| , [11|, [12|] of 
Gribov-copies effect, in which the "fc-bc" comparison was used to assess the Gribov ambiguity. A 
detailed check of Zwanziger's conjecture on the weakening of the Gribov problem with an increase 
of the lattice volume (using the SA vs. OR "one-copy" comparison) requires further studies both 
for smaller and larger lattices. 

In Fig. ^ a scatter plot is shown of the ghost dressing function for a broad range of momentum, 
combining data obtained for 7 configurations on a 56^ lattice, 14 configurations on a 64^^ lattice, 3 
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Typical snapshot of ghost propagator 

L = 56, SA vs OR on Sthconf 
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Figure 1: Ghost propagator results: the influence of the gauge fixing algorithm on the propagator, calculated 
for a typical single gauge field configuration (left), on the averaged propagator (right). 



configurations on a 72"* and 3 configurations on a 80^ lattice, all thermalized at j8 = 5.7. The gauge 
field configurations were produced with a heat-bath algorithm applying 0(1000) thermalization 
sweeps in between. We consider only momenta surviving a cylinder cut with = 1 [16]. For 
all lattice sizes GF has been carried out with the SA algorithm. Surprisingly, all the values for the 
ghost propagator fall perfectly on one universal curve (within 1% accuracy), besides those for the 
2 smallest momenta. The results, especially those found at L = 80, show that a true IR exponent K 
cannot yet be defined or does not exist at all. This is at variance with the asymptotic DSE prediction 
K = 0.595 1^ and also with K = 0.2 motivated by thermodynamic considerations in JTs]]. The 
latter estimate of K is based on the required cancellation of gluon and ghost contributions to the 
pressure, that otherwise were building up a Stefan-Boltzmann law, in the confinement phase. Note 
that downward deviations of the data at lowest momentum for each physical box size V from the 
infinite-volume curve of G{q^) are predicted by the DSE approach on a finite torus [^. However, 
we do not find such deviations. 



4. Gluon propagator results 

Fig. ^shows data for the gluon propagator computed for three different lattice sizes at j8 = 5.7. 
There, the gauge was fixed with the SA algorithm. At the present stage, the data favor a non- 
vanishing gluon propagator at zero momentum, as there is no sight of a different behavior even 
at the largest lattice volume available to us. Also, the decrease of the zero-momentum propagator 
D(0) upon increasing the volume seems to become less with bigger V . The 64^ data for D{q^) 
resemble the pattern of overshooting deviations from an universal function of momentum known 
from the DSE solutions on a finite torus though. If the lowest two or three momenta were 
removed from the plot, the picture would be less convincing in favor of an emerging plateau. Better 
statistics and data on even larger symmetric lattices (with L > 80) will make us more confident in 
this. 
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Ghost dressing function 
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Figure 2: Scatter plot for the ghost dressing function from different lattice sizes and configurations, gener- 
ated at j3 =5.7 and gauge-fixed with SA. 
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Figure 3: The gluon propagator from different lattice sizes at j3 = 5.7. The data points drawn at = 0.001 
represent the zero-momentum gluon propagator D{Q). 



5. Discussion 

We conclude that using the SA technique for the purpose of gauge fixing considerably facil- 
itates simulations of ghost and gluon propagators in Landau gauge on large lattices. We find that 
SA-based computations of the ghost propagator seem to be less affected by statistical fluctuations 
compared to other calculations where OR is used. Additionally, estimates of the ghost propagator 
at low momenta are systematically lower than those obtained after simple OR. 

A continuous decrease in slope in the ghost dressing function below 0.4 GeV does not conform 
to a simple power-law ansatz. Therefore, any attempts to extract infrared exponents from lattice 
data seem to be premature at the present stage. Qualitatively, the same behavior is seen for the 
case of SU (2) (see [^]) and also in the DSE solutions on a torus [^. However, the effects of finite 
volumes are much less than expected from there though. The same we find for the gluon propagator 
which we cannot confirm to be infrared-decreasing even at volumes larger than (13fm)^. 
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Future analytical (DSE and renormalization group) studies and lattice simulations, option- 
ally including also Z(3) flip operations into the GF procedure [0], at even larger volumes will 
hopefully help to resolve or explain the existing discrepancies between the lattice and analytical 
findings. 
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